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Abstract 

We study the adiabatic limit in the density matrix approach for a quantum 
system coupled to a weakly dissipative medium. The energy spectrum of the quan- 
tum model is supposed to be non-degenerate. In the absence of dissipation, the 
geometric phases for periodic Hamiltonians obtained previously by M.V. Berry are 
recovered in the present approach. We determine the necessary condition satisfied 
by the coefficients of the linear expansion of the non-unitary part of the Liouvillian 
in order to the imaginary phases acquired by the elements of the density matrix, 
due to dissipative effects, be geometric. The results derived are model- independent. 
We apply them to spin | model coupled to reservoir at thermodynamic equilibrium. 
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1 Introduction 



Since the fascinating work by Berry in '84|]TJ showing the existence of geometric phases 
(path-dependent phases) in vector states driven by adiabatic periodic Hamiltonians, au- 
thors in the literature have looked for geometric phases in other physical contexts. In 
particular, Joye et al. || and Berry|| independently showed that the transition proba- 
bility of instantaneous eigenstates of non-real Hamiltonians in the non-adiabatic regime 
gets an imaginary phase. This imaginary phase was measured by Zwanziger et al. for a 
two- level system 0]. The appearance of an imaginary correction to the geometric phase in 
quantum models coupled to dissipative media has been discussed in the literature [J5], ||, |7fl . 
In those references, the non-unitary evolution of the quantum system is implemented by a 
phenomenological non-hermitian Hamiltonian. This phenomenological approach has been 
extensively applied to the study of the properties of open quantum systems || |9] . 

Recently we have considered a spin | model in the presence of an external magnetic 
field and coupled to a weakly dissipative medium this system precesses with constant 
angular velocity around a fixed axis. We applied two Lindblad operators to represent the 
non-unitary part of the Liouvillians of the quantum system in contact with two distinct 
reservoirs. Using the master equations for these models, we concluded that the geometric 
phases acquired by the spin ~ instantaneous eigenstates of the Hamiltonian were not 
modified by the presence of the dissipation. The effective result of the interaction of the 
quantum system with the reservoir is the shrinking of the Bloch vector, which can be used 



to give a geometric description to the density matrix [IT], [12|, [L3|, |14 |. 

Certainly it is still an open question whether our results in reference |TJJ| are of gen- 
eral nature or particular to the dissipative models studied. We remind that they are in 
opposition to the ones derived by the non-hermitian Hamiltonian approach |5], |6|, [?|. In 
order to prove that the results of JlJJ are valid in general for reservoirs at equilibrium, we 
must have a model-independent approach. We should say that master equations are the 
proper way to study an open quantum system, whereas phenomenological non-hermitian 
Hamiltonians are supposed to give a "bone fide" description of an open system only when 
the coherence of states is not destroyed by the interaction of the system with its neighbor- 
hood. We point out that the general expression for the imaginary phase acquired by the 
non-hermitian Hamiltonians is given by eqs. (75) and (76) of reference J5|, from which we 
conclude that this imaginary phase has be geometric for any non-hermitian Hamiltonian, 
as derived in references ||, and 0] for some specific models. However, by the end 
of reference || Garrison and Wright affirm that their result should be checked out by a 
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density matrix approach. 

In reference [TO] the density matrix already appears as our central object. However, 
even there our explanation is based on the time evolution of the instantaneous eigenstates 
of Hamiltonian. In the present paper, we rederive the evolution of a quantum system 
in the adiabatic approximation directly in the density matrix formulation, which is the 
natural approach in the study of quantum systems coupled to dissipative environment. 

The Adiabatic Theorem discussed in references |15|, [16|, [l7j applies to quantum systems 
driven by unitary evolution. It states that if initially the system is in an eigenstate of the 
Hamiltonian, its time evolution is fasten, at each time, to an instantaneous eigenstate of 
Hamiltonian with the same original quantum numbers. Therefore, it is a natural choice to 
write our density matrix in the basis of the instantaneous eigenstates of the Hamiltonian 
and derive its adiabatic limit. 

The dynamics of the operator density p(t) is given by a Liouville-von Neumann equa- 
tion. In the Liouvillian we add a non-unitary term to take into account the interaction 
with a dissipative medium, that is 



-^-l = -i[H(t),p(t)]+C D p(t), (1) 

where H(t) is the time-dependent Hamiltonian of the quantum system and Co is a super- 
operator that acts on p(t) and is responsible for the non-unitary evolution of the quantum 
system. 

In the basis of the instantaneous eigenstates of the Hamiltonian H (t) (i.e., H (t) \uf, t) = 
Ei(t)\uf, t)), the dynamical equations of the elements of the density matrix are 



- i(Ei(t) — Ej{t))pfj{t) + (ui]t\£ D p(t)\uj]t), (2) 

where p^(t) = (ui,t\p(t)\uj;t). Not all elements of matrix pfj{t) are independent, due to 
the constraints that Tr[p(i)]= 1 and pfj(t) be a hermitian matrix. 

The time-dependent Hamiltonian that drives the quantum system takes into account 
its effective interaction with its neighborhood. This interaction is realized through a set 
of time-dependent external parameters R(t) = (R^(t), R 2 (t), - ■ • , R m (t)). By hypothesis, 
the time-dependence of Hamiltonian comes only from the external parameters (H(t) = 
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H(R(t))). For a periodic Hamiltonian, there are in general two distinct time scales: the 
time scale associated to transitions between two instantaneous eigenstates i and j 
of the Hamiltonian (being typically of the same order of magnitude) and the time scale 
(period) T, associated to the periodicity of the external parameters (R(T) = R(0)). 

Our aim is to derive the dynamical equations satisfied by the elements pfj(t) in the 
limit Tij <C T — > oo [|l^, |Tj], [17|] . We consider a quantum system with non-degenerate 



energy spectrum. To present the main details of the calculation, in section 2 we take the 
simplest situation when the system is not coupled to a dissipative medium, in which case 
we have to recover the known results of literature J^j. In section 3, this non-degenerate 
quantum system is coupled to a weakly dissipative medium and the density matrix has 
to be used to describe the quantum behavior of the system. In subsection 3.1, we discuss 
the conditions under which the quantum system could acquire an imaginary geometric 
phase due to the presence of dissipation. In section 4 we apply the results of section 3 to 
the spin | model coupled to two particular reservoirs at thermodynamic equilibrium to 
verify the nature (time-dependent or path-dependent) of their imaginary phases. Finally, 
in section 5, we present our conclusions. 



2 The Quantum Systems in the Absence of 
Dissipation 

In the absence of dissipation, the time evolution of the quantum system in the adiabatic 
limit is completely described by its vector state. The adiabatic evolution of vector states 
is well done in references |15| , |I~6| , [T7| . The appearance of the real geometric phase in the 
unitary evolution of vector states is beautifully described in reference |J. 

Our aim in this section is to consider a simpler physical situation to present the main 
details of the calculations of getting the adiabatic limit directly from the matrix density 
approach. This have the advantage of checking the correctness of our results. 

Eq. (|[) in the absence of a coupling with a dissipative medium becomes 



dt 



i\{Ei(t) - Ej(t)) +i(-(ui-,t 



d\uj]t) 
dt 



+ (Uj\ t 



d\uj] t) . 
dt 



^ EM-EM Plj {) + ^ Pil {t} - 



EjW-Etit) 



(3) 



Since we intend to get the adiabatic limit of the elements of the density matrix, we 
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proceed as usual [1], [L^, and factorize the dynamical phase in p^(t), 



(4) 



We distinguish three types of dynamical equations for the elements Pijit): 
i) for the diagonal elements pa(t), where i — 1, 2, • - • , N — 1, namely, 



dpu{t) 
dt 



i-1 



= E 

AT 

+ E 



1 



1 



2Re 



(^;ti[^]k i ;t) e i/ o dt ' (Ei(t ' ) -^ (i ' )) ^(*)) a 



+ 



i=i+l 



E^-E^t) 



2Re 



^;t|[^]|n,t)e-^o^(^(0-^(0)^ (t) j. (5a) 



for the elements Pij(t), where i = 1,2, ■ ■ ■ , N — 1 and j = % + 1, • • • ,N — 1, namely, 



dpijjt) 
dt 



i ^ d\uj,t) .Auj-t) 
\Ui,t\ j4 _ + {Uf,t\ 



dt 



dt 



Pijit) + 



+ 



+ 



+ 



V ^ Ui '^ dH M ) \\ u ^ t ) -% /„* dt' (SKt')-Si(t')) = . m 



+ 



Z=j+1 
i-l 



r dH{t) 
[ dt 

Ei(t) - Ei(t) 



+ 



V ^ *1 [ 1 \ u i'i f ) r i dt'mW-Eitt')) ,~ (f) y 



+ 



(5b) 



ra) for the elements Pijv(t), where i = 1, 2, • • • , JV — 1, namely, 
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dp iN {t) 
dt 



+ 



+ 



+ 



..d\ui)t) , d\u N ;t) 
{Ui]t\ — — {u N ;t\ 



kCTlM ; /« dt > ( Ei (t>)-E N (t>)) ~ (f) 

^ Ek(t)-E N (t) Pll[t) 



dt 

\dH(t) ] 



dt 



PiN{t) + 



+ 



N-l 



V ^ t ^ dH dt ) ^ t ) r -i f* dt'(E l (t')-Edt')) ~ at (f) 

4f m-m 



+ 



i-1 



v-A ; 1 1 [ d } ] \u N -t) f t dt , {Ei {t , ) _ En (t , ) ) , 
^ E N {t) -E.it) {Ph[t)) + 



N-l 

+ £ 



E N (t)-E t (t) 

\dH{t) 



(Uj',t\[ M ]\uNy^_ e ifodt'( i E i (t')-E N (1/)) 

Ei(t)-E N (t) 



(5c) 



Eqs. (5a) and ( pq) can be written as Volterra integral equations [19| of the first type, 
whereas eq. ([5c]) can be written as a Volterra integral equation of the second type. 

Following references [15] and |18|] we redefine the time scale using the period T of the 
external parameters 



(6) 



In this new variable s, the integral equations obtained from eqs. (|5]) are: 
i) for the diagonal elements f>a(s), where i = 1, 2, • • • , N — 1, 



Pii{S) 



i-1 



A' 



Pn{0) + 
1 



=i+l 



+ y / 

^ /„ K,(V 



1 



-2Re 



-2Re 



ii) for the elements Pij{s), where i = 1, 2, • • • ,N — 1 and j = i + 1, • • • , iV — 1, 
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Pij(s) = p tj (0) + 



+ 



Pij(s') + 



o i=1 



/ t\rdH(s')-\\ /\ 



+ 



+ / ^'E 



O i 



r/ * ( 7b) 



iii) for the elements pny(s), where i = 1, 2, • • • , N — 1, 



Pi7v(0) 



ds Kttijs I — — {u N ;s 



ds> 



ds' 



PiN{S 



/,... yir ^O iL.,,. _/\ 



z= 

AT-1 



i-1 



+ / ^e ^^fe^ eiTJ °' ' fa "( g '('")-^('"))(p K ( a /) ) * + 



Z=l 



Ends') -Eds' 



+ / ^'E 



<^|[^]| w '> 



iT/ Q s d S "(E N ( S »)-E 1 («"))p. ! ( s '') + 



e iT/ 8 'd5"(S i (s")-^ JV (s")) 



^(s')-^(s') 



(7c) 



Some of the integrals on the r.h.s. of eqs. (^) are of the following type 
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Iu{s) = ds 



t (^i? S '\ } ] l^j S> ) AT f,f d3"(Ei(8")-E,(8")) n( c'\ 



Ei(s') - Eds 



with i ^ Z and p(s') representing one element of the density matrix. 
By defining 



(8) 



(s)= f'ds'iEiW-Ets')), 
Jo 

we recognize Iu(s) as the Stieltjes integral PJ 



(9) 



Iu(s)= ds'Fu(s')(F 9a ^ so that Ft 



0*2 ( s ') 



(10) 



and g'jz(s) = ^(f) (we are not summing over the indices i and I). 

Assuming that Fu(s') is a piece-wise continuous function in the interval [0, s], then the 
Riemann-Lebesgue Theorem [^TJ gives the value of Ia(s) in the adiabatic limit (T — > oo) 



lim / ^(sV^'^ifs') = 



(11) 



which means that for large but finite values of T, the r.h.s. of eq. fllOD can be written in 
terms of the inverse powers of T. Integrating eq. (|I0|) by parts, we get 



hi(s) = ^ 
1 

It 



Lfti(s) 



^(g) c tT 9 ,,,( S ) 



ls=0 



_ _ / ds V T * 
iT 



il(sf) d(F u (s')/g u (s')) 



ds' 



=0 11 



(12) 



From the result ( |T2"D we obtain that the elements p>ij(s) can be expanded in powers of 



oo 

pn( s ) = Yl ( 



(13) 



n=0 



where p^ys) is the coefficient of order (^) n in the expansion. 



8 



We are interested in the adiabatic limit (T — ► oo). This correspond to substitute ex- 
pansion ( |13D in eqs. (0) and keeping only terms of order (^)° in the differential equations. 
At this order, the equations become: 

i) for the diagonal elements pu(s), i = 1, 2, • • ■ , N — 1. 

^ = =► (${*) = f$(0). (14) 

This result gives us the meaning of the Adiabatic Theorem in the density matrix approach: 
the population of an instantaneous eigenstate of the Hamiltonian does not change in an 
adiabatic process. 

ii) for the elements p~ij(s), where i — 1, 2, • • • ,N — 1 and j = i + 1, • • ■ , N. 



ds 

whose solution in the variable t is 



-K.|^ + <^|^>;- ( ,; (15) 



p*{t) = J r { " '-^-J e ->J^t'mt')-EAt')) p H {Q y (16) 

Global phases do not contribute to the density matrix. Eq. ( |16|) is compatible with eq. 
(1^) since the element pfj{t) is only different from zero if in the initial state we already 
have a superposition of the i th and j th instantaneous eigenstates of H(t). Finally, the 
difference of phases, either the geometrical and the dynamical phases, appears naturally 
in this formalism. 

We notice that eq. (|16|) includes the elements pf N {t), i = 1, 2, ■ • • , iV — 1, as well 
as the solution (13) for the diagonal terms. The inhomogeneous term in eq. ( ffc[) gives 



a contribution of order (i); this is why it does not contribute to the dynamics of the 
elements pjjy(s) in the adiabatic limit. 

The argument of the first exponential on the r.h.s. of eq. (^) is equal to the difference 
of the geometric phases acquired by the instantaneous eigenstates i and j of the Hamilto- 
nian H(t). For any closed path in the R-space, each of those phases depends only on the 
path followed by the R-parameters. As should be, in the density matrix formalism we 
recover the known results derived for the unitary evolution of instantaneous eigenstates 
of Hamiltonian in the adiabatic limit liTTl . 
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3 The Quantum System in the Presence of Weak 
Dissipation 



In general we are interested in studying a quantum system that is part of a whole system 
whose sub-systems interact with one another. This interaction allows the sub-systems to 
have exchanges among themselves. The traditional way to study a part of the whole sys- 
tem is taking a partial trace over all degrees of freedom of the complementary sub-system. 
These complementary degrees of freedom are called environment. In this approach we have 
an effective Hamiltonian that drives the dynamical evolution of the quantum system under 
study and at the same time the non-unitary part of the Liouvillian takes into account its 
interaction with the environment. In the general case the effective Hamiltonian depends 
on a set of time-dependent classical parameters. 

In reference |L0| we questioned whether the presence of dissipation could introduce an 
imaginary geometric phase in systems driven by periodic Hamiltonians. We considered a 
two-level model in the presence of two distinct Lindblad operators representing reservoirs, 
and concluded that in those cases the imaginary phases are not geometric and, based on 
eq. (19) of reference ||10|| , we affirmed that the nature of the imaginary phase depends 
on its origin. Certainly, this last affirmation has to be confirmed by a model-independent 
approach. 

The last term on the r.h.s. of eq. (§) introduces the effects of dissipation in the 
dynamics of the quantum system. We assume that the dissipation is weakly coupled 
to the quantum system. In general, the non-unitary part of the Liouvillian is written 
in a time-independent basis. Let {\vi)} be this time-independent basis and, for weakly 
dissipative interaction, we have [f22|, [23f 



where pi m (t) = (vi\p(t)\v m ) . The coefficients c^(t) take into account the characteristics 
of the environment (in the case of a quantum system at thermodynamic equilibrium 
with a reservoir composed of an infinity set of harmonic oscillators, the coefficients c^(i) 
take into account the distribution of frequencies, etc, but are time-independent). In the 
most general case of coupling with a dissipative medium we can be time dependent and 
consequently those coefficients can vary in time. 

From the beginning we chose to write the density matrix in the basis of the instanta- 
neous eigenstates of Hamiltonian H(t). In this basis, the non- unitary part of Liouvillian 




(17) 
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is written as 



■t\C D p(t)\ Uj ;t) = J2 €l(t)P?m(t)- 



The coefficients c] J m (t) on the r.h.s. of eq. (|18D are obtained from c^it) by making a 
similarity transformation in each of its indices, that is 



4L(t) = (wi;*l«ii)<wi;*l«i 3 )^3u(*)<«i a |wi;*X«i4|wm;*), ( 19 ) 

where implicit sum over the indices l\, I2, h and U is meant. We point out that even in the 
case when the coefficients cj^ are time independent, their analogous in the instantaneous 
basis can acquire a time dependence through the transformation (jl9|). From the scalar 
products in the transformation ([Li]) the coefficients cf m (t) do not get a dependence on the 
variation of Hamiltonian or of any other external classical parameter. 

Following the same steps as we did in section 2, we obtain the dynamical equations 
for the elements p~ij(s) (see their definition in eq. @) in the limit T — > 00: 

i) for the diagonal elements pu(s), % — 1, 2, • • • , N — 1. 



dpu(t) 
dt 



N-l 



(<#(*(*)) - 4^(Mt)))p«W + 4* (*(*))■ (20) 



In writing the r.h.s. of eq. ( P0{ ) we have already implemented the condition Tr(p(t)) = 1. 
In the most general case the dissipation couples the dynamics of the diagonal elements 
of the density matrix. From eq. (|20|) we recover the adiabatic behavior of the quantum 
system in the presence of the dissipative medium when the non-unitary part of the Li- 
ouvillian has null diagonal terms in the instantaneous basis (cjj =0). In this situation 
the population of an instantaneous eigenstate of Hamiltonian does not vary along the 
adiabatic process. In this case the quantum system does not transfer to the environment 
energy due to electronic transitions. 

The elements of matrix C(t) are defined as 

c u (t) = 4(t)-4 N (t), (21) 
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i, I — 1, 2, • • • , N — 1. The general solution of eq. ( p0|) is 



T ( e fZdt'C(t') 



dt' 



T ( e J°dt"C(t") 



J jk 



(0+^(o) 



(22) 



The operator T means a time-ordering operator [21]. From eq. (|1^) we obtain that 
the coefficients c\l(t), i,l = 1, 2, ■ ■ • , N, are the same for any basis of the instantaneous 
eigenstates of Hamiltonian. This fact avoids any ambiguity in the imaginary phases in 
the time-ordering terms. In the general case the imaginary phases in eq. fl2"2] ) are time- 
dependent. We postpone the discussion under what conditions the time-ordering integrals 
in eq. fl2"2] ) can be rewritten as a path integrals in a suitable parameter space. 

From eq. (^) we recover solution ([TJP in the absence of dissipation (c|| = 0, i, I — 
1,...,N). 



ii) for the elements Pij(s) 



1,2,-.. ,N 



1 and j — % + 1, 



N. 



dpijjt) 
dt 



dt 



+ 



(23) 



{l,m} 



where the set of pair of indices {/, m} are those which satisfy: Ei(t)—E m (t) = Ei(t)—Ej(t). 
To simplify our discussion, we order the instantaneous eigenenergies such that: if Ei(t) < 
Ej(t) then i < j. Since the indices (i, j) of the elements of the density matrix in eq. (|23|) 
are chosen such that i < j, we get that the elements pi m (t) that contribute to the r.h.s. of 
this equation are such that / < m. Once the energy spectrum of the quantum system is 
non-degenerate, if the elements pi imi (t) and p~i 2 m 2 (t) that contribute to r.h.s. of eq. (|23|) 
are distinct then we must necessarily have l\ ^ I2 and mi 7^ m 2 . 

Let us suppose that the dynamics of M elements f>ij(t) are coupled by the presence 
of dissipation and their dynamics are given by eq. (|23|). Due to the fact that each pair 
(/, m) is unique we may relabel them by using only one index: (Z i5 rrii), i — 1, 2, • • • , M. 
Eq. (p3|) is rewritten as: 



dt 



{Uk;t\[-^\u k ;t) ) + (u mt ;t 



M 



\U,. 



Pi, 



At) + ^2^l k Pl k m k (t), 



k=l 



(24) 
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i = 1, 2, M. We distinguish two possible situations: 



1) the dynamics of the elements Pij{t) are not coupled by the coupling to a weakly dissi- 
pative medium: M = 1. 



In this case, eq. (|24|) reduces to 



dt 



(ui-, t\ (^k; *)) + ( u j; t\ (^K; *)) + CjiPijit)' ( 25 ) 



whose solution is (see eq. 



(t) = e -ift<#(EiW-Estt')) e -fi<# 



Pg(0). (26) 



Here again the coefficients c^(i) that appear in the last phase on the r.h.s. of eq. (|26"D are 
independent of the chosen basis of the instantaneous eigenstates of Hamiltonian. In the 
general case this phase is time-dependent. In next sub-section we discuss the conditions 
satisfied by the coefficients c^(t) in order to this integral becomes path-dependent. 

Eq. ( p6|) reduces to eq. QlGf) in the absence of interaction with a dissipative medium. 

2) the dynamics of M elements p%j{t) are coupled due to the presence of dissipation: 
M > 1. 

Making the change of variables 



PkmM = e 



i , d \ u l ■ !*') , ,, d|u m . 



(27) 



where z = 1, 2, • • • , M, eq. fl24[) becomes 



dt 



A ik (t)p hmk {t). 



The elements v4jfc(t) of matrix A(t) are defined as 



(28) 



Wfc'*! — dt 7 — Wi-M — dp — {v>m, k ,t\ 2T nu mv t\ ^7 



(29) 



The solution of eq. ( ^8f) is 
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L \ / J tj 



T being the time-ordering operator |24j . 
From eqs. (01) and (p7|) we finally have 



(30) 



K;f| 



d|n mi ;t') 



P?m,-(0)- 



X 



(31) 



By choosing a new basis of eigenstates of Hamiltonian the elements A^it) get an irrelevant 
real phase that does not contribute to the average value of any physical operator. As in 
the two previous discussion, the time-ordering integral that appears on the r.h.s. of eq. 



( j3ll) due to the presence of dissipation is time- dependent. In the next sub-section we give 
the conditions necessary for this imaginary phase to be path-dependent of a suitable set 
of time dependent parameters. 

Eq. flilf) is recovered from eq. fl31"|) in the absence of a coupling with a dissipative 
medium. 



3.1 Conditions to Obtain an Imaginary Geometric Phase 



In the most general case, eqs. (|2"2"D, ( |2"6| ) and fl3"T| ) do not give an imaginary geometric 
phase correction to the real Berry phase [I]. In this subsection we discuss the mathematical 
requirements that the coefficients c^ k (k(t)) have to satisfy in order to those imaginary 
phases are geometric (path-dependent). 

Differently from the real phases acquired by the evolution of the instantaneous eigen- 
states of Hamiltonian in the adiabatic limit [|l^, 16, [L7[], the imaginary phases in eqs. fl22[ ), 



(Pq) and (31) have no ambiguity due to the arbitrariness of the basis of the instantaneous 
eigenstates of Hamiltonian. Consequently in the case of the imaginary phases acquired 
by the elements of the density matrix p H (t) due to the presence of dissipation, they can 
be written as an integral over a suitable time- dependent parameter space if we have an 
integration over t of a function f(t) that has the general form 



/(t) = X>(t) 



dt 



(32) 
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with the functions ^(i) satisfying the following conditions: 

1) the functions *&i(t) are not explicitly time-dependent; 

2) the time-dependence of functions ^(t) come only from their dependence on the set of 
parameters k(£) = (ki(t), k2(t), ■ ■ ■ , hit)). 

We point out that we do no restrict the regime of the time variation of the set of 
parameters k(t) and it has not to be a periodic function in time. 

We begin by discussing the imaginary phase on the r.h.s. of eq. (p2[). On the r.h.s. of 
eq. (^) we have a time-ordering integrals of a matrix C(t) whose elements are (see eq. 

my- 



C u (t) = 4(t) - c% N (t), i, l = 1, 2, • • • , N - 1. (33) 

The relation between c\\(t) and the coefficients of the non- unitary part of Liouvillian 
written in a time-independent basis is given by eq. ([TOD 



c u(t) = (uh^i^iu^tlv^cll^ivi^^ivijuut}, (34) 

where i, I = 1,2, ■ ■ ■ , N. In order to be able to write the time-ordering integrals on the 
r.h.s. of eq. (|22|) as a path-dependent integral the coefficients cj^(i) must have the form 



and the functions ^i s ijf) have to satisfy the two previous conditions mentioned. We 
stress out that the path in the k-parameter space has not to be closed. 

Due to the hermicity property of the operator Cop{t) the elements of matrix C (see 
eq.(|2"ID) are real. If condition fl3"5] ) is satisfied the geometric phase introduced by the 



presence of dissipation is purely imaginary. 

Once eq. (^) is valid for the coefficients c^Jt), the non-diagonal elements pf m {f) 
have two distinct situations in the limit T — >• oo: 

i) for M = 1, the elements pf m {t) get an imaginary geometric phase. From the hermicity 
property of operator C D p(t) we obtain that c*^(t) = ^c^(t)J which leaves open the 
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possibility that the coefficient c^(t) has an imaginary part. For master equations where 
cMt) has an imaginary part, the dissipation gives a real correction to the Berry's phase]!]]. 
This correction to the Berry's phase has not to be a closed-loop in the k-parameter space. 

ii) for M > 1, the elements pf^it) get an imaginary geometric phase due to the interaction 
with a weakly dissipative medium. The presence of the dissipation couples the dynamics 
of different elements of the density matrix and the integral over the k-path involves a 
matrix that in general does not commute with itself at different instants. This k-path has 
not to be a closed loop. As in case i if the coefficients c^(i) have an imaginary part they 
give a correction to the real geometric phase [H. 

From eq. ( p5| ) we see that the presence of the dissipative media gives an imaginary 
geometric phase only if the coupling between the quantum system and its environment is 
externally driven by a varying external parameter. In both situations the time-variation 
of the set of the parameters k(t) has not to be slow. 



4 Adiabatic Limit of the Density Matriz of 
the Spin | Model Coupled to a Reservoir 

For systems in contact with reservoirs at thermal equilibrium whose the weak coupling 
constant does not vary in time, the effect of the presence of the dissipation is to destroy 
the coherence in the quantum system in a time-dependent process. That is the case of 
two distinct reservoirs coupled to a spin | model studied in detail in reference |L0 . 



By using directly eq. fl35|) we want to verify in the next sub-sections if the conclusions 



about the nature of the imaginary phases in reference [10] are corrected or not. 



4.1 Dephasing Process in a Two-Level System 



An interesting process well studied in the standard textbooks |23|, is the phase destroy- 
ing process which might appear due to elastic collisions. We consider a spin 1 variable 
(two level model) coupled to a time dependent magnetic field precessing around the z-axis 
with uj constant precession frequency. The external magnetic field B(t) has norm B and 
makes a 6 angle with the z-axis. 

For the sake of later calculations, it is convenient to define two unitary transformations: 
the first one, R(u,t), takes us to the rotating frame where the Hamiltonian is no longer 
time dependent; the second one, D(B,8,uj), diagonalizes the effective Hamiltonian (time 
independent) that drives the dynamics of the final matrix representation of the density 
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operator. We call this the diagonal frame ]TIJ . 

The master equation to this process written in the diagonal frame is 



d k 

-j t PD if) = -i [\ x a z , p D (t)] + - (cr z p D (t)a z - p D (t)) 



(36) 



where k is the dissipation constant at zero temperature. The weak coupling regime is 
characterized by the condition ^ <C 1 with Ai = yj (pB cos(#) — ^ ) 2 + p?B 2 sm 2 (8). 

We define p H (t) to be the density matrix in a basis of the instantaneous eigenvectors 
of Hamiltonian. The relation between p (t) and po{t) is Hj 



p H (t) = Vi(t)T>p D (t)T>V(t), 



where the matrix V(i) is equal to 



V(t) 



|)e 2 _ sm (j) e 
sm(|Je 2 cos^je 



and 



2 
2 



(37a) 



(37b) 



D = —=Vi + A ct z + -^=Vl - A cr a 
y/2 V2 



(37c) 



with A = (ypB cos(#) — ^\ and ±Aj are the eigenvalues of the effective Hamiltonian. 



The master equation (|36|) in the instantaneous basis of the Hamiltonian for arbitrary 
value of us is 



d_ 
Jt' 



p H (t) = -i (pB + ^)a z -^-a n (t),p H (t) +^-C%p H {t) 



to 



(38a) 



where 



C H D p H {t) = A 2 a n (t)p H (t)a n (t) + KVl^[e- wt (a n {t)p H {t)a + {t) + a + (t)p H '(t)cr n (t)) + 
+ e^(a n (t)p H (t)a4t) + <r_(t)p H (tK(t))] + (1 - A 2 ) [e- 2i ^a + (t)p H (t)a+(t) + 
+ e 2i "V_(t)p ff (t)a_(t) + a + (t)p H \t)o_(t) + <x_(t)p H (i)a + (t)l - p H (t). (38b) 
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In eqs. (pH), we use the definitions: 



a n (t) = 



cos(#) 
— sin(6>)e 



- cos(0) 



(39a) 



<T+(t) 



sin(fl) 



cos 2 (f)e" 



■iuit 



sm 



sin(6>) 



(39b) 



and <J-{t) = \ o+{t) 

From eqs. (^) and (p6| ) we see that to obtain the nature of the imaginary phase (if they 
are time or path-dependent) we need the coefficients c^J, and cj^ gotten from C^p H (t) in 
the adiabatic limit (u — *■ 0). However, the weakness condition on the dissipation constant 
k (j^) does not impose any constraint to the ratio K The interesting case happens when 
the adiabatic evolution and the dissipation effect are of the same order: - ~ 1. 



In the adiabatic limit, weak dissipation regime and for - ~ 1, eq. (|38b| ) give us 



CO 



and 



and 



ii 



cs = k O(-) 



(40a) 



,,12 
'21 



(40b) 



We neglect the contribution of terms of order kO(-^) to eqs. (^) and (26) since they 
are of the same magnitude as the terms of higher order in (i). The constant k can not 



be written as eq. (pq) and consequently the imaginary phase due to the weak coupling 
to the dissipative medium is not geometric. For completeness we substitute the values of 
the coefficients in eqs. ( p2|) and (pB]), and obtain [p6[ 



(41a) 



e -2i)iBt e ~iuj(l-cos(8))t e ~kt „H 



pf 2 (o). 



(41b) 



From eq. ( |41a| ) we get that the process continues to be adiabatic even in the presence 
of dissipation. On the r.h.s. of eq. (|41b|) the first phase gives the difference of the 
dynamical phases of the eigenstates of Hamiltonian, the second phase gives the difference 
of the geometrical phases acquired by the instantaneous eigenstates of H(t) and the last 
phase is the time-dependent phase due to the coupling of the quantum system with the 
dissipative medium. Its effect is to destroy the quantum coherence in the system. 
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4.2 Adiabatic Limit of a Two-Level Model in Thermal 
Equilibrium 

In our next application, we consider the spin ~ model under the influence of an external 
magnetic field as described in section 4.1 but now coupled to a reservoir of electromagnetic 



fields at thermal equilibrium [10 . The master equation of this model in thermodynamic 



equilibrium in the diagonal frame is [23, 25] 



— p D it) = -i [AiCTg, p D (t)] +k(n + l) (2cr_p jD (t) a+ - p D (t) cr + a_ - ct + (J-Pd if)) + 
+kn {2a + p D (t) ct_ - p D (t) a^a + - a^a + p D (t)) , (42) 

where k is the dissipation constant at zero temperature and n is the average number of 
excitations of the weakly coupled thermal oscillators at inverse temperature j3. 

By doing the transformation (|3"7|), the master equation ( pE2"D in the instantaneous basis 
of the Hamiltonian for arbitrary value of uj becomes 



dt 



p H (t) 



pB 



a z -^a n (t),p H (t)j + k£%p H (t) 



(43) 



where 

c H DP H (t) 



2n + l 



2p H (t) + (1 - A 2 ) o n (t)p H {t)a n (t) 



-2iuit 



o + (t)p H (t)o + {t) 



2iut 



a4t)p H (t)a.(t) +(1 + A 2 ) a + (t)p H (t)a„(t) + a4t)p H (t)a + (t) 



- AVT^A* [e iut (a n (t)p H (t)a_(t) + o_(t)p H (t)o n (t)} + 



+ e 
1 
2 



^[a+(t)p H (t)a n (t) + a n (t)p H (t)a + (t)^ 
{p H (t),Aa n (t) + v / I^A 2 (e-^a + (t) + e*"V_(t)} + 2A a + {t)p H {t)a^{t) 



o_(t)p H (t)a + (t)\ + VT^[e-^(a n (t)p H (t)a + (t) - o + (t)p H (t)a n (t)) + 



+ e^(a4t)p H (t)a n (t) - a n (t)p H \t)a.(t) 



(44) 



a n (t) and er+(t) are given by eqs. ( |39a| ) and ( |39b| ) respectively. 

As in sub-section 4.1, we consider the adiabatic limit and the weak dissipative regime. 
In those regimes we consider the case - ~ 1 and the coefficients c}^ and cj^ gotten from 
C^p H (t) (see eq. ([£§)) are: 
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2(1 + n) +0{ 



UJ 



and 



c 22 



2n + 0{ 



(45a) 



and 



-.12 
'21 



k - (1 + 2n) + . (45b) 

By the reason discussed in sub-section 4.1 the terms k C(^) do not contribute to the 
adiabatic limit. For the present model we see that neither the coefficients cj^ and cj^ 
can be written as eq. ( j35|) and consequently the phases due to the coupling of the spin | 
model with a reservoir of electromagnetic fields at thermal equilibrium are not geometric. 
The solutions of eqs. ( j2~2"|) and (|26|) for this model are [^6] 



PnW= Pn(0) 



n 



l + 2n 



e -2k(l+2n)t _j_ 



T7 



l + 2n 



(46a) 



and 



p g(t) = e -^Bt e -Ml-cos(e))t e ~k(l + 2n)t p H^ y ^ 



From eq. (|46aj) we get that due to the weak interaction with the dissipative medium, we 



have a very slow phenomena but the adiabatic character of the quantum process is lost. 
As in the previous sub-section, the constants k and n can not be written as eq. (^) and 
consequently the imaginary phase due to the weak coupling to the dissipative medium 
is not geometric. From eq. (|46b ) we obtain that the interference effects between the 



instantaneous eigenstates of Hamiltonian is destroyed by the coupling with the dissipative 
medium. This effect depends on the dissipative constant k and the time elapsed. 

5 Conclusions 

We discuss the behavior of the density matrix of non-degenerate quantum systems whose 
dynamics are driven by a master equation in which the unitary part is a periodic Hamil- 
tonian with period T (T —>■ oo). Our discussion is model- independent. We recover the 
known results of the literature pD] for the quantum systems that are not coupled to any 
external dissipative medium. The difference of geometric phases necessary for measuring a 
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physical effect associated to the existence of these phases appears naturally in the density 
matrix approach. The obtainment of the adiabatic limit for the density matrix in a closed 
system is used as a simple situation where we present the details of the calculations. 

Next we consider the quantum system coupled to a general weak dissipative medium. 
We take the case - ~ 1 when the effects on the dynamics due to the slow evolution 
and the dissipative attenuation are of the same order. From eq. (^) we obtain that the 
adiabatic nature of the phenomena disappears if the non-unitary part of the Liouvillian 
has non-null diagonal terms in the instantaneous basis of Hamiltonian. 

In the general case the loss of coherence due to the coupling to a weak dissipative 
medium is a time-dependent phenomena (see eqs. ( p6|) and (fH])). In eq. (^) we have 
our main result where we obtain the condition that the coefficients cj^ (t) have to satisfy 
in order to the imaginary phase be geometric. This only happens when the interaction 
between the quantum system and its environment is time-dependent and is a consequence 
of the time variation of a set of parameters. As in the case of imaginary geometric phases 
appeared in the transition probabilities of non-real Hamiltonians in the non-adiabatic 
regime 0, [| here also the integral over the parameter space k does not have to be a 
close loop. The condition ( |35| ) allows us to say directly from the expression of C% if the 
dissipative effect in the master equation gives an imaginary geometric correction to the 
Berry phase Q. In the case where condition (|3f|) is fulfilled and the coefficients cf m (t) has 
an imaginary phase, the presence of dissipation gives a correction to the real geometric 
phase. 



The condition (|35D is not satisfied by the coefficients c^^(t) in the linear expansion of 
the non- unitary part of the Liouvillian (see eq. ([18])) that represents the weak interaction 
of a non-degenerate quantum system with a reservoir at thermodynamic equilibrium. This 
implies that for those type of couplings the imaginary phases are time-dependent. In the 
present work we extend the validity of the results derived for two particular dissipation 
mechanisms discussed in reference [|T0j that are rediscussed in section 4 using the approach 
presented in this work. 

The non-hermitian parts in the Hamiltonians in references |5|, [], [jj take into account 
the losses of a quantum system to its environment (a suitable reservoir of degrees of 
freedom at equilibrium). In these models, the couplings between quantum systems and 
their neighborhoods in the ll ab initio" Hamiltonians do not depend on the time variation 
of a set of classical parameters. This fact, togheter with result ( P5| ) (the condition to the 
existence of an imaginary geometric phase) put in check the correctness of the imaginary 
geometric phases in the literature due to dissipative effects H H, EH- m or der to verify if 
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the imaginary phase for a quantum model described by a phenomenological non-hermitian 
hamiltonian truly exists — being of true geometric origin, and not a fake one due to eqs. 
(75) and (76) — one must apply the approach presented here. The presence of imaginary 
phases of geometric and time-dependent natures in a quantum system yield distinct rates 
of loss of coherence, and such distinction is experimentally detectable. 

Finally we point out that we have applied the results derived in here to quantum 
systems in contact with reservoirs at thermodynamic equilibrium, but they apply equally 
well to models whose interaction with the environment varies in time. Only in this 
situation the imaginary phase acquired by loss of coherence can be eventually geometric 
(if and only if the condition (^) is fulfilled). 
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